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ABSTRACT 

The  influence  of  sound-speed  fluctuations  on  propagation  of 
a CW  signal  in  an  ocean  with  a uniformly-sloping  bottom  and 
a horizontal  surface  is  analyzed  using  ray  theory.  The  mean 
sound-speed  structure  is  modeled  as  bilinear,  with  bottomed  source 
and  receiver  above  and  below  the  SOFAR  axis,  respectively.  The 
horizontally-independent  fluctuations  oscillate  with  a 12-h 
period  in  the  upper  ocean.  An  examination  is  made  of  possible 
types  of  rays  for  down-slope  propagation  that  might  exist, 
depending  on  bottom-slope  angle  and  source-receiver  separation. 

The  total  acoustic  field  is  investigated  for  its  dependence  on 
these  parameters  and  time.  For  certain  conditions  when  up  to 
three  rays  comprise  the  mean  total  field,  three  patterns  of  time 
evolution  are  described,  each  of  which  may  have  significant 
amplitude  variations.  Numerically-computed  examples  of  each  type 
are  presented.  The  linear  relationships  between  phase  variations 
of  individual  rays  and  the  sound-speed  fluctuations  are  derived. 
Then,  formulas  are  developed  to  explain  the  most  frequent  behavior 
of  the  relative  amplitude  and  phase  of  the  multipath  total  field. 
Predictions  from  the  formulas  show  very  good  agreement  with  the 
numerical  calculations. 


INTRODUCTION 


This  paper  continues  an  investigation,  begun  in  Ref.  1,  into 
effects  of  certain  combinations  of  ocean-environmental  features  on 
CW  acoustic  propagation.  Reference  1 was  concerned  with  the 
influence  of  a combination  of  spatially-uniform  currents  and  upper- 
ocean,  tidal  sound-speed  fluctuations  on  deep-ocean  transmissions 
over  a horizontal  bottom.  Under  the  model  assumptions,  it  was 
found  that  total-field  phase  is  affected  by  both  current  and 
sound-speed  fluctuations,  and  that  the  latter  more  significantly 
contribute  to  amplitude  variations.  The  small  number  of  other 
acoustic  studies  in  which  combinations  of  oceanic  variations  have 
been  considered  are  discussed  in  Ref.  1. 

The  primary  purpose  of  this  paper  is  to  examine  effects  of 

sound-speed  fluctuations  on  down-slope  acoustic  propagation  over 

a nonhorizontal  bottom.  One  model  of  interest  is  a uniformly- 

sloping  bottom,  which  may  be  regarded  as  grossly  approximating 

portions  of  the  topography,  at  some  locations,  from  continental 

or  island  shelves  to  the  abyssal  ocean.  For  simplicity  we  focus 

on  this  model,  rather  than  possibly  more  generally-applicable 

ones  consisting  of  several  piecewise-linear  segments.  We  note 

that  our  procedures  permit  extensions  to  such  cases.  Another 

reason  for  considering  the  uniform  slope  is  that  bottom 

variation  has  been  excluded  from  most  previous  published 

propagation  studies.  Exceptions  concern  analyses  of  properties 

2 

of  refracted/bottom-reflected  rays  and  surface-ref lected/bottom- 

3 4 

reflected  rays  over  a uniform  slope.  Another  investigation 

showed  how  to  analyze  refracted/bottom-reflected  rays  propagating 

in  an  arbitrary  direction  over  a uniform  slope,  by  relating  such 
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problems  to  ones  in  an  ocean  with  a uniform  horizontal  sound- 
speed  gradient.  Environmental  fluctuations  were  not 
considered  in  these  studies.  Moreover,  their  principal 
objective  was  relatively  near-shore  propagation,  for  which  a 
linear  or  constant  mean  sound  speed  is  appropriate.  In  the 
present  case,  for  transmissions  to  depths  below  the  SOFAR 
axis  (at  about  1000  m) , a bilinear  mean  sound-speed  profile  is 
physically  more  reasonable,  as  well  as  analytically  tractable 
(see,  for  example.  Refs.  5,  6,  and  7).  In  particular,  we 
assume  that  both  source  and  receiver  are  bottomed,  above  and 
below  the  SOFAR  axis,  respectively.  In  this  paper  we  study 
acoustic  propagation  along  a range  normal  to  the  constant-depth 
lines.  This  situation  is  a good  approximation  of  some 
experimental  conditions,  is  illustrative  of  sloping-bottom  effects, 
and  is  relatively  easy  to  analyze.  Other  propagation  directions 
could  be  studied  by  applying  and  extending  the  procedures  in 
Ref.  4.  It  is  also  assumed  that  for  given  bottom-slope  angle  a, 
the  source-receiver  distance  R is  short  enough  so  the  receiver 
depth  does  not  exceed  the  abyssal  ocean  depth  (about  4700  m) . 

In  order  to  understand  the  characteristics  of  the  chosen  model, 
we  investigate  intervals  of  values  for  a (between  7 and  15  degrees) 
and  R (5  and  20  km),  which  include  those  typical  of  relatively 
short-range  transmissions  from  off-shore  sources. 

In  Sec.  I the  sound-speed  fluctuations  are  modeled  by 
pivoting  sinusoidally,  with  semidiurnal  period,  the  upper-ocean 
sound-speed  gradient  about  the  surface  sound  speed.  The  range  R 
is  assumed  short  enough  to  ignore  horizontal  variations. 

Motivation  for  this  model  in  terms  of  the  semidiurnal  tide  is 


3 


discussed  in  Ref.  1 and  references  therein.  We  remark  that  the 
model  usefulness  is  not  necessarily  restricted  to  this 
particular  form  of  fluctuations.  In  Sec.  II  the  wide  variety  of 
possible  types  of  RBR  and  SRBR  rays  that  may  occur  are 
considered.  It  is  found  that  the  possibilities  in  the  sloping- 
bottom  geometry  change  dramatically,  depending  on  values  of 
parameters  including  a and  R.  Conditions  are  discussed  so  that 
the  received  total  field  consists  of  no  more  than  three 
particular  rays,  and  such  cases  are  treated  subsequently. 
Expressions  are  given  for  the  amplitude  and  phase  of  these  three 
rays  at  the  receiver. 

Numerical  calculations  for  the  relative  phase  and  amplitude 
of  the  total  field,  as  functions  of  time,  are  discussed  in 
Sec.  III.  Three  distinct  patterns  are  found  for  the  multipath 
cases,  each  of  which  may  have  significant  amplitude  variations. 
All  show  the  interaction  of  the  fluctuations  and  the  sloping 
bottom,  and  for  one  type  in  particular,  a consequence  is  the 
disappearance  and  later  reappearance  of  a ray.  In  Sec.  IV 
relations  between  the  phase  perturbations  of  individual  rays  and 
the  sound-speed  fluctuations  are  obtained.  Then,  formulas  are 
derived  for  total-field  variations  in  terms  of  ray-phase 
perturbations,  for  the  most  commonly-occurring  type  of  behavior 
of  the  multipath  total  field.  Comparisons  between  predictions 
from  the  formulas  and  computations  in  Sec.  Ill  show  quite  close 
agreement,  in  cases  where  conditions  hold  for  validity  of  the 
formulas.  Main  results  of  the  paper  are  summarized  in  Sec.  V. 
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I.  MODEL  FORMULATION 

We  model  a variable-depth,  off-shore  oceanic  region  by 
assuming  that  the  ocean  surface  is  flat  and  that  the  ocean  bottom 
is  a sloping  plane  intersecting  the  surface  at  an  angle  a. 

The  line  of  intersection  is  regarded  as  the  coast  line  of  our 
model.  We  use  a Cartesian  coordinate  system,  with  origin  at 
some  point  on  the  coast.  Increasing  values  of  x and  z 
indicate  distance  normal  to  the  coast  toward  the  deep  ocean  and 
oceanic  depth,  respectively,  and  the  ocean  surface  is  the  plane 
z = 0.  The  y-direction  is  parallel  to  the  coast  so  as  to  produce 
a right-handed  coordinate  system.  The  ocean  bottom  in  this  system 
is  the  plane  2 = x tan  a.  As  discussed  in  the  Introduction,  we 
restrict  attention  in  this  paper  to  propagation  normal  to  the 
coast,  thereby  eliminating  dependence  of  both  the  range  and  the 
sound  speed  on  the  coordinate  y.  The  sound  source  is  placed  on 
the  ocean  bottom,  at  depth  dQ,  and  the  receiver  is  also  on  the 
ocean  bottom,  at  depth  d^.  The  horizontal  distance  between  source 
and  receiver  is  the  range  R and  may  be  expressed  as 

R = (d  - d )cot  a . (1) 

*\  o 

Figure  1 illustrates  the  geometry  of  the  propagation  range  and 
the  parameters  introduced  above,  along  with  some  others  to  be 
discussed  later.  We  note  that  our  solution  procedures  could  be 
adapted  to  analyze  other  hydrophone  locations  as  well  as 
propagation-range  orientations.  One  additional  constraint  on 
our  particular  model  is  that  the  sloping  bottom  should  not 
extend  to  a depth  d0  of  the  abyssal  ocean  beyond  the  land-mass 
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rise,  i.e.  d_  < d . Otherwise,  a piecewise-linear  bottom 

K B 

model,  in  which  the  slope  at  depth  d becomes  parallel  to  the 
ocean  surface,  may  be  more  appropriate. 

We  choose  for  our  model  of  sound  speed  c the  expression 

c(x,z,t)  = Cg ( 1+ f ( z ) ] + c(x,z,t),  (2) 

where  cg  is  the  fixed  sound  speed  at  the  ocean  surface  above 
the  source,  i.e.  at  x = dQ  cot  a,  and  t is  time.  The  function 
f(z)  describes  the  depth  variation  of  sound  speed  in  the  absence 
of  6,  which  incorporates  changes  due,  for  example,  to  tidal 

effects.  The  dimensionless  quantities  f and  c c"1  satisfy 

- -1  . . 
c c <<  f <<  1. 

We  specify  the  form  of  the  mean  variation  f(z)  as  a 
bilinear  model,  of  a type  used  previously  (for  a recent  example, 
see  Ref.  1) , to  reasonably  approximate  measured  deep-ocean 
profiles  while  simplifying  the  acoustical  analysis.  We  thus 
select 

f ( z ) = -zgjCg1  » 0 £ z £ dA  (3a) 

and 

f (z)  = [cA  - Cg  + g2(z~dA^CS1  ' dA  - z - X tan  01  * 

In  Eqs.  (3) , cA  is  the  value  of  the  mean  sound  speed  which 
occurs  at  the  SOFAR  axis  depth  d^.  The  quantities  g-^  and  g2  are 
the  magnitudes  of  the  sound-speed  gradients  above  and  below  the 
SOFAR  axis,  respectively.  We  may  now  express  the  mean  sound 
speed  at  the  sound  source  c0  and  at  the  receiver  cR  as 


- — — 


Cs  - Vo 


CR  ’ CA  + WV  ' (4b» 

where  we  have  chosen  to  consider  the  case  d < d,  < d„ 

OAR 

exclusively  in  this  paper. 

Typical  parameter  values  appropriate  for  studying  sound 

O 

propagation  initiated  near  Bahaman  shores  of  the  Sargasso  Sea 
are  c_  = 1541  m sec  c.  = 1487  m sec  * and  d„  = 1000  m, 
leading  to  the  value  0.0540  for  g , and  g^  = 0.0146.  For 
definiteness,  we  will  use  these  fixed  numerical  values  in 
computations  throughout  this  paper.  Nonetheless,  we  seek  to 
construct  a model  with  more  generality  and  wider  possible 
application  than  just  to  this  oceanic  region  or  to  any  specific 
acoustic-propagation  range.  Thus,  we  will  present  results  for 
source-receiver  distances  R between  5 and  20  km,  and  for  bottom- 
slope  angles  a in  the  range  between  7 and  15  degrees.  These 
values  include  those  typical  of  relatively  short-range 
transmissions  from  off-shore  sources.  in  this  paper  we  will  only 
discuss  investigations  for  a fixed  source  depth  dQ  = 300  m. 
Results  for  this  case  were  regarded  as  typical  of  source  depths 
not  near  either  the  surface  or  SOFAR  axis.  For  this  fixed  value 
of  d_  and  given  a and  R values,  it  follows  that  d„  is  deter- 
mined  by  Eq.  (1) . For  a value  of  d = 4700  m,  requiring  the 
receiver  remain  on  the  slope  and  below  the  SOFAR  axis  imposes 


700  m < R tan  a < 4400  m 


(5) 


I 
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We  note  that  other  fixed  parameter  values  and  intervals  of  R and 
a values  can  be  examined  by  further  application  of  our  solution 
procedures. 

We  consider  next  the  effect  of  tidal  variations  on  sound 
speed  via  c.  An  earlier  study*  considered  the  tidal  effects  on 
sound  speed  in  the  deep  ocean,  using  observations  of  temperature 
fluctuations  (as  in  Ref.  9),  an  analysis  of  state  equations  to 
relate  sound-speed  fluctuations  to  those  of  temperature,  and 
the  assumption  of  sufficiently  short  propagation  ranges.  The 
deep-ocean  model  in  Ref.  1 includes  an  oscillation  with  period 

A 

12  h of  both  the  depth  dft(t)  of  minimum  sound  speed  and  the 

A 

gradient  above  dA<t)  , with  the  surface  sound  speed  remaining 
fixed.  The  semidiurnal  component  is  a strong  feature  of  temperature 
fluctuations  nearer  to  shore  as  well,  and  frequency  spectra  of 
currents  there  indicate  even  more  concentration  of  energy  in 
this  component*®  Since  we  are  aware  of  no  evidence  to 
invalidate  the  deep-ocean  model  for  c in  the  continental-slope 
regions,  we  shall  employ  it  here.  We  note  that  in  any  application 
of  our  model,  more  general  time  variation  could  be  incorporated 
if  observations  warranted.  Moreover,  certain  other  forms  of 
vertical  structures,  either  singly  or  in  combination,  could  be 
treated  by  the  method  in  Sec.  IV. 

Therefore,  we  set 


c(z,t)  = 


Uzsin  (irt/6)  , 0 < z < d.  (t) 

— — A 

A 

0 , z > dA(t) 


where  p is  the  magnitude  of  the  perturbation  gradient,  measured 
in  sec  *, 


and  t is  measured  in  hours  so  that  for  convenience  the 
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tidal  perturbation  at  t = 0 is  absent.  The  time-dependent  SOFAR 
axis  depth  is  closely  approximated  by1 

dA(t)  * d^fl+y  (91+92}  lsin  • (7a) 

/>  A 

It  follows  that  the  time-dependent  sound  speed  c (t)  at  d (t)  is 

A A 

A A 

cA(t)  - c(dA(t),t)  = cA  + Asin  (irt/6)  , (7b) 

where  the  non-negative  parameter 

A = ydAg2(g1+g2)  1 (7c) 

/v 

is  the  amplitude  of  the  sound-speed  fluctuation  of  d (t)  . Values 
of  A up  to  0.30  m sec  1 will  be  examined  here,  as  in  Ref.  1. 
Equations  (6)  and  (7)  imply  that  the  maximum  variation  of  sound 
speed  with  time  occurs  at  the  mean  SOFAR  depth  d and  is  equal  to 

A 

Ag2  (g1+g2) . Figure  2 illustrates  the  model  by  showing 
representative  sound-speed  structures  corresponding  to  maximum, 
zero,  and  minimum  perturbations. 


II.  RAY  ANALYSIS 


We  assume  that  acoustic  propagation  occurs  under  conditions 
for  which  ray  theory  is  valid.  In  this  section,  we  discuss  the 
ray  geometry  and  amplitude  and  phase  of  any  ray  in  the  variable- 
depth  channel  of  Sec.  I with  no  sound-speed  perturbations.  Thus, 
the  sound  speed  is  given  by  Eq.  (2)  with  c = 0 and  Eq.  (3),  and 
parameter  values  are  as  in  Sec.  I.  It  follows  that  ray  consist 
of  piecewise  circular  arcs.  Moreover,  because  the  vertical  plane 
containing  the  source-receiver  direction  is  perpendicular  to  the 


- 
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coast  line  of  our  model,  rays  remain  entirely  in  this  plane; 

otherwise,  rays  making  bottom  reflections  would  be  deflected 

4 

out  of  this  plane.  If  9 is  the  inclination  angle  of  a ray 
measured  positive  clockwise  from  the  horizontal,  then  the 
usual  form  of  Snell's  law, 

c cos  0 = constant,  (8) 

is  valid,  although  the  constant  changes  after  each  bottom 
reflection. 

We  now  consider  determination  of  the  rays.  In  order  to 
obtain  analytical  ray  descriptions  which  are  also  numerically 
accurate,  we  begin  by  hypothesizing  a possible  ray  path.  For 
our  model,  rays  can  be  characterized  by  specifying  a pair  of 
integers  [K,L] , where  K and  L are  the  number  of  surface  and 
bottom  reflections,  respectively.  We  note  that  under  our 
particular  model  assumptions,  it  can  be  shown  by  geometrical 
considerations  and  numerical  calculations  that  there  exists  at 
most  one  ray  for  each  pair  of  integers  [K,L] . Next,  the  equations 
governing  the  ray  are  determined,  which  in  general  form  a 
nonlinear  system  for  the  ray  angle  at  the  source,  0q,  and  other 
parameters.  Numerical  root-finding  procedures  are  then  used  to 
determine  values  for  the  variables.  The  validity  of  each  solution 
is  established  by  testing  that  the  ray  does  not  pass  outside 
the  ocean  on  its  way  to  the  receiver. 

The  simplest  conceivable  ray  has  K - L = 0.  Its  governing 
equations  are  Eqs.  (1) , (4) , (8) , and 

R = gI1(cAtan  eA  - cQtan  0q)  + g^c^an  0ft  - cRtan  0R)  , (9a) 


•whmmmm 
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where  0 and  0 are  (positive)  ray  angles  at  the  SOFAR  axis  and 
A R 


receiver,  respectively.  Equation  (9a)  is  derived  in,  for  example, 


Ref.  4.  The  [0,0]  ray  is  illustrated  in  Fig.  1.  The  angles 


0,0,  and  0 are  shown  and  are  all  positive.  The  requirement 
O A R 


that  this  ray  never  reaches  the  surface  may  be  expressed  as 


0 > - cos  1 (c  c ■*■)  . 

o os 


Additional  inequalities  which  are  imposed  for  given  a and  R 


are  Eq.  (5)  and 


0R  - “ ' 


which  precludes  acoustic  energy  from  reaching  the  receiver  through 


the  ocean  bottom. 


Numerical  results  for  the  existence  of  the  [0,0]  ray  are 


shown  in  Fig.  3.  Values  for  0Q  may  be  found  by  iterating  for 


tan  0 using  Newton's  method  for  most  (ot,R)  pairs  in  the  Figure. 


Equation  (10)  is  satisfied  only  above  and  to  the  left  of  the 


curve  labeled  3 in  Fig.  3.  Moreover,  Eq.  (9b)  is  violated  below 


curve  8.  Since  the  first  inequality  of  Eq.  (5)  is  not  satisfied 


in  the  shaded  area  above  curve  1,  it  follows  that  the  [0,0]  ray 


exists  only  in  the  region  of  Fig.  3 bounded  by  curves  1,  3,  and  8, 


We  remark  here  that  the  second  inequality  of  Eq.  (5)  fails  in  the 


other  shaded  area,  below  curve  7. 


Another  ray  of  interest  has  K = 1 and  L = 0,  and  is  also 


illustrated  in  Fig.  1.  The  equations3  determining  the  geometry 


of  the  [1,0]  ray  are  Eqs.  (1),  (4),  (8),  (10), 


" " --  ■ ■ ” 


u ,,  M .IT***:,-, 
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R = gT1  (c.tan  0.  - c tan  0 - 2c_  tan  0_) 

1 A A o os  s 

+ g”1  (c^tan  0 - c tan  0 ) , 


(11a) 


and 

0 < - cos_1(c  c”1)  , (lib) 

O O o 

where  0g  is  the  (positive)  ray  angle  at  the  surface.  As  illustrated 

in  Fig.  3,  the  [1,0]  ray  is  found  by  numerical  calculations  to 

exist  in  the  region  of  the  ja,R)  plane  bounded  by  the  curves 

labeled  2,  4,  and  8.  Below  curve  4,  the  [1,0]  ray  fails  to 

exist  because  Eq.  (10)  is  violated,  while  below  curve  8,  the 

[1,0]  ray  disappears  because  Eq.  (lib)  is  violated. 

The  third  type  of  ray  of  principal  interest  in  this  paper 

is  that  with  K = 2 and  L *=  1.  As  illustrated  in  Fig.  1,  its 

bottom  reflection  is  at  depth  dR^,  where  the  bottom  incident  ray 

angle  is  0 . and  the  sound  speed  c is,  from  Eqs.  (2)  and  (3a) , 
ol  B 1 

CB1  * CS  ' VB1  • (12a> 

Additional  conditions  needed  to  specify  the  [2,1]  ray  are  Eqs.  (1), 
(4),  (8),  (10), 

R * g.1(c  tan  0 -2c_(tan  0_.  + tan  0oO)  - c tan  0 
1 A AS  SI  S2  o o 

+ cB1ltan  001  - tan (2a-00l) ] } + g~1(cAtan  0A~cRtan  0R) , (12b) 

(dBfdo)cota  * g^1(cBitanGBr2cstanesrcotan0o)  * (12c) 

and 

2a  - 0B1  < " cos’1(cBlCS1)  ‘ 


( 12d) 
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The  angles  0gl  and  0g2  are  the  (positive)  ray  angles  at  the  first 


and  second  surface  reflections,  as  shown  in  Fig.  1.  Note  that 
Snell's  law  provides  two  independent  equations,  since  the  bottom 
reflection  changes  the  constant  in  Eq.  (8). 

The  numerical  calculations  for  the  [2,1]  ray  were  performed 
by  first  estimating  tan  0^  and  c^.  Then,  Snell's  law  was 


used  to  determine  0^  and  0g^,  after  which  Newton's  method  was 


applied  to  Eqs.  (12a)  and  (12c)  to  converge  to  the  correct  value 


for  cBl  for  the  estimated  value  of  tan  0fil.  Values  were  obtained 


for  0_-,0-,  and  0_  from  Snellte  law,  and  tan  0_.  was  re-estimated 
S2  A R B1 


from  Eq.  (12b) . This  process  was  continued  until  convergence 

was  obtained.  Efficient  calculations  for  this  ray,  as  well  as 
for  others,  relied  on  accurate  initial  estimates.  It  was  found 
that  the  [2,1]  ray  exists  over  most  of  the  (a,R)  plane  in  Fig.  3, 
although  not  below  curve  5 because  Eq.  (10)  is  violated.  In 
addition,  we  assumed  that  the  source  hydrophone  is  directional, 
transmitting  significant  sound  energy  only  for  launch  angles  0Q 
satisfying  |0Q|  less  than  some  value  0qD-  This  can  be  accomplished, 
for  example,  by  using  an  omnidirectional  "point"  source  with  a 
reflector,  as  described  for  example  in  Ref.  11,  or  an  array  of 


omnidirectional  sources.  With  the  value  0qD  = 40  , the  [2,1] 


ray  is  restricted  to  the  left  of  curve  6 in  Fig.  3. 

For  certain  parameter  values  considered  here,  rays  of  type 
[0,L]  exist,  with  L ^ 1 bottom  reflections  which  must  lie  above 
the  SOFAR  axis.  As  indicated  in  Fig.  3,  we  ha  e shown  that 
various  members  of  this  family  of  rays  occur  for  the  shallow 
receiver  positions  above  curve  2.  When  these  rays  exist,  they 
could  in  principle  be  significant  contributors  to  the  total  received 
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field,  since  they  suffer  relatively  small  spreading  and  bottom 
loss  because  of  their  short  path  length  and  near-grazing  bottom 
incident  angles.  On  the  other  hand,  for  the  latter  reason, 
such  rays  are  particularly  sensitive  to  actual  nonlinear  bottom 
topography  in  any  specific  experiment.  For  our  investigations 
of  total  field  in  later  sections,  we  choose  parameter  values  for 
which  rays  of  this  type  do  not  occur. 

We  remark  that  the  combination  of  bilinear  sound  speed  and 
sloping  bottom  in  our  model  permits  a remarkably  rich  variety  of 
rays.  Figure  3 suggests  the  patterns  that  can  arise,  but  depend- 
ing on  assumptions  and  parameter  values,  the  situation  can  become 
more  complicated.  For  instance,  if  0oD  is  larger  than  40°,  not 
only  does  the  [2,1]  ray  exist  over  more  of  the  (a,R)  plane  shown, 
but  also  new  rays  may  occur.  The  [3,2]  ray  may  be  present, 
since  its  initial  angle  ranges  from  0q  = - 40.5°  for  (a,R)  = 

(1° ,11  km)  to  0Q  = ~ 85.0°  for  (a,R)  = (15°, 5 km),  for  the  (a,R) 
values  in  Fig.  3.  Indeed,  other  rays  with  more  than  two  surface 
reflections  could  also  occur.  The  importance  of  these  steep  rays 
to  the  total  field  would  depend  critically  on  such  factors  as 
bottom  loss  and  scattering.  As  another  example,  a piecewise-linear 
bottom  dramatically  changes  the  ray  possibilities;  more  than  one 
ray  sometimes  occur  for  given  numbers  [K,L]  of  surface  and  bottom 
reflections. 

In  view  of  the  foregoing  considerations,  for  subsequent 
analysis  of  the  total  received  field,  we  shall  restrict  (a,R) 
values  to  the  regions  in  Fig.  3 where  some  or  all  of  the  [0,0], 
[1,0],  and  [2,1]  rays,  but  no  others,  exist.  We  henceforth  denote 
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these  as  the  n = 1,  2,  and  3 rays,  respectively.  Note  that  for 
the  rays  considered  in  this  paper,  we  can  identify  n with  K + 1. 

The  travel  times  Tn  for  rays  n = 1 and  n = 2 can  be  adapted 
from  Ref.  5,  while  that  for  ray  n = 3 can  be  calculated  using 
procedures  of  Ref.  5.  The  results  are 

T.  = g ,3ln[c  (l+sin0_ ) c * (l+sin0  ) 3]  + 

11  o A A o 

g^ln  [cR  (l+sin0A)  c”1  (l+sin0R)  _1]  , (13a) 

T2  = g”1lnIcQ(l-sin0s) (l+sin0A)cA1(l+sin0o)  1 

x (l+sin0c)-1]  + g3  ln[c  (l+sin0  )c  ^(l+sine  ) 3]  , (13b) 

b fc  K A A K 

and 

T3  = g^1  ln{cQ (l-sin0gl) (l+sin0Bl) (l-sin0g2) (l+sin0A) 

x c"1 (l+sin0Q) -1 (l+sin0sl) _1 (l+sin0s2) _1 [l+sin(2a-0B1) 1 _1} 

+ g^1  In [cR(l+sin0A) c”1 (l+sin0R) -1]  . (13c) 

These  expressions  can  be  evaluated,  using  Eqs.  (1) , (4) , (8) , (12a) , 
and  (12c)  where  appropriate,  once  the  source  angles  0Q  of  the 
three  rays  have  been  determined  numerically. 

The  formulas  for  the  geometrical  spreading  loss  can  be 
developed  from  results  and  methods  in  Ref.  2: 

A.  = Rc2  sin0_  sec30  I tan0  IgT1  (c^cote.  - c-1cot0  ) 

X O R O*  O 1 A A O O 

+ g21(CAlcot0A  " CRlcot8R)]| 


9 


(14a) 


15 


A_  = Rc2  sin0_  sec30  tan0  [g, * (c. ^cotO  - 2c  1 cot0c 
2 o R o o 1 A A S S 

- c 1 cot0  ) + g *(c  * cot©.  - c 1 cot©  ) 1 , (14b) 

o o 2 A A R R 

and 

A3  = R sin0R  sec0o(g'1{c21sec20Bltan0o(c-Jcot0B1-c;1  cot0Q 

- 2c'1  cotOg)  + [cot(2a-0B1)  - cot0R1]  (dcBl/deo) } 

+ c0lsec2 (2a-0B1>  t (l+tan(2a-0Bl)cot0Bl)  (dcfil/deo) 

- cn1tan0  tan(2a-0o1 )cot0  ] IgT1 [c -1cot0 

Bl  o B1  B1  I 1 A A 

- ZCg1  cot0s2  - Cg^cot (2a-0Bl) I + g21(cA1cot0A 

- cRX  + (1-R^lT1)  cos0B1sec(0Bl^-2a)J  . (14c) 

In  Eq.  (14c)  the  derivative  term  is 

dc  ,/d0  = c tan0  secO  (csc0  ,-csc0  -2csc0ol ) (cote., -cota)  \ 

Bl  O O O O Bl  O SI  Bl 

( 14d) 

and  the  horizontal  distance  between  source  and  bottom  reflection  is 

**81  = ^i1(C0-CBl)COt°  * (14e) 

We  assume  the  sound  source  emits  a harmonic  signal  with 

circular  frequency  w and  with  maximum  amplitude  unity,  i.e. 

-2 

1 dyn  cm  at  distance  1 m.  It  follows  that  ray  n arrives  at  the 

receiver  in  the  form  A„  sin  (iot  - 4>  ) . The  amplitude  A is 

n n n 

affected  by  source  directivity  and  by  the  physical  processes  of 
spreading  loss,  boundary  loss,  attenuation,  and  scattering,  but 
the  latter  two  effects  are  excluded  here.  At  any  surface 
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12 

reflection,  the  model  of  Mackenzie  leads  to  no  amplitude  loss 

and  a phase  shift  of  it.  The  loss  and  phase  shift  <J>3  at  the 

bottom  reflection  of  ray  3 may  also  be  computed  from  formulas 

in  Ref.  12.  Two  parameters  characterizing  bottom  properties 

appear  in  the  formulas  for  B3  and  4> 3 . The  first  (called  a/0  in 

Ref.  12)  is  the  ratio  of  complex  to  real  parts  of  the  bottom 

attenuation  factor.  Although  this  quantity  in  general  depends 

on  acoustic  frequency  and  bottom  composition,  a frequency- 

12 

independent  value  may  be  derived,  and  we  choose  a value  of 
0.0158  for  a bottom  which  has  been  used  for  the  Sargasso  Sea 
area.^  The  second  parameter,  o,  is  the  impedance  ratio  of  the 
bottom, 


o 


(pxcx) (PB1CB1) 


(15a) 


where  the  corresponding  bottom  density  p and  sound  speed  c are 
Px  = 2.00  gm  cm-3  and  c^  = 1757  m sec  3 . An  approximate  formula 
for  water  density  PB1(in  gm  cm-3)  at  depth  dR1(in  m) , which  can 
be  derived  from  Eq.  (20)  of  Ref.  14,  is 


PB1  = 1.0243  + (8.2716  x 10"6)dB1  . 


Thus,  the  amplitude  of  ray  n can  be  written  as 


A 

n 


\~1/2  B D ( 0 ) , 

n no' 


(15b) 


(16a) 


if  we  define  B^  = B^  = 1.  The  directivity  factor  D is  normalized 
to  a maximum  value  of  1.  We  note  that  our  definition  of  6 „ 
requires  D(0  ) to  be  significantly  less  than  1 for  |0  I > 0 

'O'  OD 

Similarly  if  <J>^  = 4>2  = 0,  the  phases  4»n  satisfy 


♦ = a)Tn  - (n-l)ir  - <f>n  . (16b) 

Using  Eqs.  (13)  and  (14),  the  values  for  the  functions  in  Eq.  (16) 

are  found  after  the  numerically-determined  quantities  such  as 

0 are  known, 
o 


III.  TOTAL-FIELD  RESULTS 

Under  the  conditions  described  in  Sec.  II  for  which  Fig.  3 

is  applicable,  the  total  field  at  the  receiver  is  given  by 

3 

A sin  (wt  - 4>)  = y)  ^ An  sin(u>t  - 4>n)  , (17a) 

n=l 

in  which  the  amplitude  A and  phase  4>  are  determined  from 
3 3 

S = ^ An  sin  <Pn  , C=  ^ An  cos  4>n  , (17b) 

n=l  n=l 

and 

2 2 1/2 

A=(S  + C)  , sin  <J>  = S/A  , cos  4>  = C/A  . (17c) 

Expressions  for  An  and  4>  are  given  in  Eqs.  (16).  The  amplitude 

A is  set  to  zero  if  ray  n is  absent, 
n 

The  preceding  section  focused  on  ray  analysis  in  the  absence 
of  the  sound-speed  perturbation  c.  Here,  we  present  results  of 
the  numerical  investigation  of  Eqs.  (17)  for  our  model  including  c. 
Under  the  assumptions  in  Sec.  I,  the  time  variation  of  the  sound 
speed  is  sufficiently  slow  that  it  can  be  ignored  during  the  time 
required  for  an  acoustic  signal  to  propagate  from  source  to 
receiver.  From  Eq.  (6)  it  follows  that  rays  in  the  perturbed 
medium  are  of  the  same  form  as  those  when  c = 0.  Moreover,  the 
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launch  angle  0 and  other  variables  for  a ray  of  type  n = 1,  2, 
o 

or  3 can  be  calculated  from  the  appropriate  governing  equations 
in  Sec.  II,  but  with  perturbed  values  of  g^,  d^,  and  c^.  Given 
values  of  a and  R,  the  diagram  for  existence  of  rays  in  the 
perturbed  medium  differs  only  slightly  from  Fig.  3,  specifically 
in  the  location  of  certain  boundary  curves.  A number  of  different 
combinations  of  values  for  parameters  a and  R were  selected,  but 
all  the  calculations  we  shall  discuss  here  used  A = 0.3  m sec 
a)  = 206  Hz,  and  D(6q)  = 1 for  |0Q|  <_  40°  and  0 for  larger  1 0Q | . 

Results  for  this  A are  presumably  typical  of  those  for  other 
(small)  values  of  A.  The  directivity  pattern  used  here  may  be 
regarded  as  an  idealization  of  more  realistic  patterns  (for 
examples  constructed  from  endfire  linear  arrays  of  isotropic 
point  sources,  see  Ref.  Is) . Time  was  varied  over  12  h since 
that  interval  is  the  period  of  c. 

From  the  many  cases  examined,  several  characteristic  types  of 
behavior  in  the  total  field  were  observed,  which  we  now  discuss. 

The  simplest  type  is  obviously  when  a and  R are  such  that  only  a 
single  ray  appears,  which  from  Fig.  3 would  be  either  the  n = 2 
or  3 ray.  Typically,  amplitude  A varies  only  slightly  with  time; 
the  maximum  variation  of  relative  amplitude  over  12  h was  found 
to  be  less  than  1 dB  in  all  single-path  cases  examined.  Phase  4> 
is  sinusoidal  in  time,  with  period  12  h.  The  peak-to-peak  phase 
variation  was  0.42  cycles  for  the  n = 2 ray  when  (a,R)  = (15°,  7 km). 

It  would  be  anticipated  that  single-path  cases  possess  a simple 
relationship  between  environmental  and  acoustical  variations. 

This  is,  of  course,  an  important  consideration  in  certain  acoustical 
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experiments  which  monitor  ocean  processes.  In  our  model,  it  is 
significant  that  single-path  cases  occur  at  relatively  longer 
ranges  than  multipath,  as  is  clear  from  Fig.  3. 

As  shown  in  Fig.  3,  other  values  of  a and  R permit  two  rays, 
specifically  the  pairs  (n=0,  n=l)  and  (n=l,  n=2) . Total-field 
amplitude  and  phase  are  illustrated,  for  an  example  of  the  former 
case,  with  (a,R)  = (11^,  7 km) , in  the  solid  curves  of  Fig.  4. 

The  amplitude,  which  is  relative  to  that  of  the  unit-amplitude 
source,  is  a smooth,  roughly  sinusoidal  function  of  time,  with 
peak-to-peak  variation  of  9.8  dB.  We  note  that  individual  ray 
amplitudes  have  very  small  amplitude  variation  with  time, 
analogous  to  the  single-path  case.  Hence,  the  conclusion  follows 
that  the  solid  amplitude  curve  in  Fig.  4 reflects  multipath  effects. 
The  phase,  shown  relative  to  1.0  cycles  for  clarity,  is  also 
smooth,  is  virtually  sinusoidal,  and  has  peak-to-peak  variation 
of  0.68  cycles.  This  pattern  was  the  most  common  one  among  our 
double-path  results,  and  we  shall  refer  to  it  as  type  I behavior 
of  the  multipath  total  field.  A theoretical  description  of  this 
case  will  be  discussed  in  Section  IV. 

Another  variety  of  double-path  behavior  is  shown  by  the  solid 
curves  for  (a,R)  = (12°,  8 km)  in  Fig.  4,  which  we  denote  by  type  II. 
The  mean  state  possesses  n = 0 and  1 rays,  but  the  n = 0 ray  is 
nearly  tangential  to  the  bottom  at  the  receiver.  The  sound-speed 
perturbation  can  cause  the  ray  to  violate  Eq.  (10) , thereby 
removing  the  ray  from  the  total  field.  This  occurs  in  Fig.  4 
from  about  t = 1.2  h to  4.8  h,  during  which  period  the  phase  and 
amplitude  display  the  characteristics  of  single-path  propagation 
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discussed  earlier.  Outside  this  time  interval,  the  double-path 
behavior  is  typical  of  type  I.  The  jump  in  amplitude  at  the 
transition  points  is  about  9.5  dB.  In  combination  with  the 
approximately  7 . 5 dB  variation  in  amplitude  from  t = 4.8  h to 
t = 9 h,  it  is  clear  that  substantial  amplitude  variation  will 
occur  in  this  example.  The  phase,  shown  relative  to  0.75,  is 
discontinuous  at  the  exit  and  entrance  points  of  the  n = 0 ray. 

We  note  that  the  actual  phase  jumps  at  those  points  are  the  differences 
shown  on  Fig.  4 (about  +0.4  cycles)  modulo  one  cycle.  Type  II 
behavior  points  out  another  effect  of  the  interaction  of  sound- 
speed  fluctuations  and  the  sloping  bottom. 

An  example  of  a third  class  of  double-path  behavior  found 
among  our  computational  results  is  illustrated  in  Fig.  5.  This 
is  an  n = 2 and  3 case  for  (a,R)  = (11°,  15  km) , with  peak-to-peak 
amplitude  and  phase  variations  of  7 . 1 dB  and  0.31  cycles, 
respectively.  However,  neither  relative  amplitude  nor  relative 
phase  of  this  type  III  case  is  close  to  the  sinusoidal  curves 
typical  in  type  I.  For  example,  amplitude  fades  occur  near 
t = 6 h and  12  h,  over  time  scaler  much  shorter  than  the  scale 
of  environmental  variation.  This  pattern  has  been  observed  in 
many  multipath  studies  (see  for  example  Refs.  1 and  16),  although 
typically  in  situations  with  more  rays  than  are  present  here. 

In  the  mean  state,  the  phases  of  the  n = 2 and  3 rays  are  nearly  ir 
radians  apart,  and  the  total-field  amplitude  is  only  about  20% 
of  the  individual  amplitudes.  Thus,  with  the  perturbation  applied, 
both  amplitude  and  phase  are  sensitive  to  small  differences  in 
the  rates  of  perturbation  of  the  phases  of  the  individual  arrivals. 
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This  sensitivity  precludes  application  of  the  theoretical  analysis 
to  be  discussed  in  Sec.  IV. 

From  Fig.  3 it  is  seen  that  for  certain  (cx,R)  combinations, 
all  three  rays  can  contribute  to  the  total  field.  Exactly  as  for 
the  double-path  calculations,  type  I behavior  was  most  frequently 
found  for  triple-path  calculations.  Relative  phase  and  amplitude 
for  a typical  case,  (a,R)  = (8°,  8 km)  shown  as  solid  curves  in 
Fig.  6,  are  very  similar  to  the  analogous  curves  in  Fig.  4. 

The  amplitude  curve  shows  a total  variation  of  7.8  dB  and  is 
somewhat  less  symmetric  about  t = 6 h than  the  corresponding 
curve  in  Fig.  4.  Peak-to-peak  phase  variation  is  0.77  cycles. 

Type  II  behavior  also  occurs  for  triple-path  cases. 

Figure  3 shows  that  either  the  n = 1 or  n = 3 rays  can  be  lost 
near  boundary  curves  3 and  6,  respectively.  We  shall  not  discuss 
cases  near  boundary  8,  where  two  rays  can  disappear.  Phase  and 
amplitude  are  shown  in  Fig.  6 for  (a,R)  = (10°,  10.5  km),  for 
which  the  n = 1 ray  is  absent  between  t = 0.5  h and  t = 5.5  h. 

In  that  interval,  the  double-path  behavior  is  characteristic  of 
type  I.  In  fact,  over  each  continuous  portion  of  the  curves, 
typical  type  I behavior  occurs.  We  remark  that  although  type-  III 
double-path  behavior  is  conceivable  when  one  ray  is  lost  from  a 
triple-path  case,  we  never  observed  this  situation.  From  Fig,  6, 
the  total-field  amplitude  variation  is  about  12.2  dB,  while  the 
changes  over  each  continuous  portion  in  this  example  are  about 
2 dB.  The  phase,  shown  in  Fig.  6 relative  to  1.1  cycles  for 
clarity,  is  approximately  sinusoidal  over  its  continuous  portions. 
As  in  Fig.  4,  the  actual  phase  jumps  could  be  larger  than 
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suggested  by  Fig.  6,  in  which  phase  modulo  one  cycle  is  shown. 
We  note  finally  that  type  III  triple-path  behavior  could 
conceivably  occur  for  three  rays  if,  for  example,  the  phase  of 
the  vector  sum  of  the  n = 2 and  3 rays  happened  to  differ  by 
nearly  ir  radians  from  the  phase  of  the  (strongest)  n = 1 ray. 
However,  this  was  not  observed  among  the  triple-path  cases  we 
examined . 


IV.  APPROXIMATE  TOTAL-FIELD  MODEL 


The  computational  results  of  Sec.  Ill  suggest  that  analytical 
approximation  of  the  total  field  phase  and  amplitude  should  be 
possible  for  behavior  of  type  I (and  for  continuous  portions  of 
type  II  cases) . The  first  step  is  to  approximately  determine 
perturbations  to  phases  of  individual  ray  arrivals  produced  by 
the  sound-speed  perturbation.  Then,  we  shall  show  how  these 
results  are  useful  in  approximating  the  total  field. 

A 

From  Eqs.  (7),  the  time-dependent  SOFAR-axis  depth  d (t) 

ft 

A 

can  be  expressed  in  terms  of  cA(t) , and  similarly  for  the  time- 
dependent  upper-layer  gradient.  Thus,  the  time  dependence  of 
the  phase  *n  of  any  ray  (n  = 1,  2,  or  3)  can  be  considered  a 
consequence  of  the  time  dependence  of  the  single  environmental 

A 

quantity  cA(t) . Define  A$n  as  the  phase  fluctuation  of  ray  n 
due  to  the  sound-speed  perturbation.  From  the  first  two  terms  of 
a Taylor  series  expansion  and  Eq.  (7b) , it  follows  that 


A*  (t)  = * (c  ) - * (c  ) 
n n A n A 


<»V8CA) 


Mt) 


(18a) 


CA  °A 


where 


A ( t)  = cA(t)  - CA  = A sin  (Tft/6) 


(18b) 


Partial  differentiation  is  used  in  Eq.  (18a)  to  indicate  that 
A<t>n  depends  on  other  variables;  just  one  term  appears  in  Eq.  (18a) 

a 

because  cA  is  regarded  as  the  only  time-dependent  environmental 
parameter.  Next,  we  use  Eq.  (16b)  to  write 


(34>  /3c  ) 
n A 


= u>  ( 3T  /3c,.) 
n A 


2 8 


n 


(18c) 


''A  A 

which  is  exact  for  n = 1 and  2 but  ignores  the  contribution 

- ( 3 ^ 3/^ CA)  | for  n = 3.  This  term  can  be  estimated,  with 
CA 

difficulty,  from  Eqs.  (15)  and  the  equations  for-  the  n = 3 ray. 
Alternatively,  its  effects  can  be  shown  to  be  negligible  by 
actual  numerical  calculations  for  the  perturbed  n = 3 ray  (for 

A 

illustration,  a typical  value  of  | (3<j>-j/3cA)  I A I is  .001  cycles 

• r* 


for  our  parameter  values  and  (a,R)  = (10°,  15  km)).  Once  an 
expression  for  the  travel  time  derivative  is  obtained,  A<l>n  is 
approximated  through  Eqs.  (18) . 

A 

Direct  calculation  of  3T  /3c  from  Eqs.  (13)  is  feasible  but 

n a 

A* 

very  unwieldy,  because  of  the  implicit  relationships  among  c and 
the  variables  describing  the  ray  geometry.  By  differentiating 
Eqs.  (13)  and  using  Eqs.  (2) , (3) , (7c) , (8) , (9a) , (11a) , (12b) , (12c) , 
and  (12d) , it  can  be  shown  after  considerable  manipulation  that 
Eqs.  (18)  yield: 

—2  —1  _i 

A-t^ft)  = sin (irt/6)  {In  [cQ (l+sin0A) cft  (l+sin0o)  ] 


- CglCj1  sin0A  - c-lsin0o) } , 


(19a) 
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A*2  (t)  « sin  ( rrt/6)  (In  [cQ  (l-sin0g)  (l+sin0A) c^ 

x (l+sin0  ) -1 (l+sin0c)  1]  - c_ (c“1sin0.-2c”1sin0 
O S S A A S 5 

- c”1sin0o) } , (19b) 

and 


A#3  (t)  = sin  (irt/6)  (ln{cQ (l-sin0gl) (l+sin0Bl) 

x (l-sineg2) (l+sin0A)c"1(l+sin0o)"1(l+sin0sl)“1 
x (l+sin0g2)  1 [l+sin(2a-0Bl>  l-1}  - cg  [c^sinO^^ 

- 2cg1sin0gl  - 2cg1sin0g2  + c^sinO^  - c^sin  (2o-0B1) 


- c"1  sin0Q] ) . (19c) 

Equations  (19)  can  be  evaluated  from  previously-computed  geometric 
quantities  of  the  unperturbed  rays. 

We  note  that  the  phase  perturbations  are  independent  of  g2 
and  0R,  parameters  characterizing  the  ray  geometry  below  the 
SOFAR  axis.  That  is,  the  sound-speed  fluctuation  produces  changes 
in  ray  geometry  below  the  SOFAR  axis,  but  these  changes  do  npt 
affect  the  phase  perturbations  in  Eqs.  (19).  To  the  order  of 
terms  retained  in  Eqs.  (19),  phase  perturbations  depend  only  on 
those  portions  of  the  rays  above  the  SOFAR  axis. 

One  way  to  present  the  results  of  Eqs.  (19)  is  to  regard 

A#  (t)  as  a function  of  a and  R and  show  the  level  curves  of  this 
n 

function.  In  view  of  Eqs.  (18) , it  is  sufficient  to  consider 
level  curves  of  8n*  Given  values  of  a and  R,  the  quantity  &n 
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could  be  obtained  by  interpolation  from  the  level  curves.  It  is 
particularly  convenient  to  plot  the  level  curves  for  n=l  and  2 
with  respect  to  R and  Rtana  = dR-dQ(by  Eq.  (1)).  To  see  this,  we 
note  that  for  any  ray,  Bn  depends  only  on  the  mean  ray  geometry 
above  the  SOFAR  axis.  Consequently,  Bn  will  be  the  same  for  all 
members  of  each  family  of  rays  with  identical  ray  lobes  for 
z < dA-  Since  the  n = 1 and  2 rays  in  particular  are  unaffected  by 
the  ocean  bottom  for  z < d^,  members  of  any  such  family  correspond 
to  rays  reaching  a receiver  at  depth  d related  to  R and  a by  Eq.  (1) . 
Because  the  locus  of  receiver  depths  is  a circular  arc,  the  level 
curves  themselves  are  circular  arcs  in  the  (dR-dQ,R)  plane.  Using 
values  for  quantities  previously  determined  numerically,  we  plot 
level  curves  for  B-^  and  62  in  this  plane  as  solid  curves  in  Fig.  7(a) 
and  7(b),  respectively.  Also  indicated  by  dashes  are  curves  1 and  2, 
of  minimum  (a  = 7°)  and  maximum  (a  = 15°)  bottom  angles,  respectively. 
In  addition,  the  curves  3 and  4,  demarcating  the  existence  of  the 
n=l  and  2 rays  as  described  in  connection  with  Fig.  3,  are  shown  in 
Figs.  7(a)  and  7(b).  We  note  that  B-^  and  B2  are  negative  for  the 
parameter  values  we  considered,  because  a net  increase  in  sound 
speed  produces  a decrease  in  travel  time  and  phase. 

The  calculation  of  B3  is  not  as  simple,  because  the  n = 3 ray 
reflects  above  the  SOFAR  axis  and  the  ray-family  characterization 
useful  for  n - 1 and  2 does  not  apply.  Level  curves  of  B3  were 
approximated  from  values  of  B3  computed  directly  from  Eqs.  (18)  and 
(19c) . The  computations  were  performed  for  values  at  the  intersection 
points  of  a rectangular  grid  in  the  (a,R)  plane.  Each  computation 
requires  determination  of  the  mean  ray  geometry,  but  the  method  is 
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straightforward.  The  results  are  shown  in  Fig.  7(c),  along  with 
the  limiting  curve  5 and  directionality  curve  6 from  Fig.  3.  The 
level  curves  have  the  same  general  shape  as  those  in  Figs.  7(a)  and 
7(b).  As  an  example,  to  compute  the  value  of  &3  for  (a,R)  = (10°, 

10  km),  we  note  d-d  = 1.76.  The  appropriate  position  is  shown 
as  a cross  on  Fig.  7(c).  By  linear  interpolation,  we  find 
63  = -1.07. 

We  now  develop  a theoretical  explanation  for  the  behavior  of 
total-field  phase  and  amplitude  for  type  I multipath  behavior  (and 
type  II  multipath  between  the  entrance  and  exit  times  of  a ray) . 

We  suppose  ♦ and  A are  the  total-field  phase  and  amplitude  (relative 
to  the  unit-amplitude  source)  in  the  absence  of  sound-speed 
perturbations.  We  define  A4>  and  AA/A  as  the  total-field  phase  and 
relative  amplitude  deviations  due  to  perturbations,  i.e. 

A 

A$  = 4’(ca)  - ♦(c  ) = 4>(t)  - $ , (20a) 

and 

AA(t)/A  = (A(cft)  - A(ca))/A(ca)  = (A ( t)  - A) /A  . (20b) 

We  note  that  AA/A  is  the  ratio  of  the  amplitude  deviation  produced 
by  the  perturbations  and  the  unperturbed  amplitude.  It  follows 
from  Eq.  (20b)  that  the  logarithmic  amplitude  anomaly  is 


AL ( t)  = 20  log1(Jll+(AA(t)/A)  ] = 20  log1QA(t)-20  log1QA  . (20c) 

Thus,  knowledge  of  AA/A  determines,  via  Eq.  (20c),  the  difference 
in  decibels  between  the  amplitudes  with  and  without  the  perturbations. 
If  AA/S  is  sufficiently  small,  then  AL  is  approximately  equal  to 
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(8.68) (AA/A) . 

We  assume  that  variations  in  individual  ray  amplitudes  from 

the  perturbations  have  negligible  effect  on  the  total  field.  This 

assumption  is  confirmed  numerically  by  results  in  Sec.  Ill,  and 

could  be  demonstrated  analytically  by  the  same  methods  used  to 

confirm  it  in  other  problems  (see,  for  example,  Ref.  17).  Thus, 

we  regard  each  Ap  as  independent  of  time,  and  we  need  only  consider 

the  influence  on  AA/A  and  A*  of  phase  perturbations  A4>  (t)  . We 

n 

define,  for  m ? n, 


e ( t)  = A4>  (t)  - A*  (t)  , 
mn  m n 


1 < m,  n < 3 , 


(21a) 


where  e is  the  phase-perturbation  difference,  expressed  in  cycles, 
between  rays  m and  n.  We  further  assume  that  the  maxima  over  t of 
any  2nem^(t),  the  perturbation  differences  expressed  in  radians,  are 
sufficiently  small  compared  to  one.  However,  even  values  of  2^0 
larger  than  0.5  give  acceptably-accurate  approximations.  Moreover, 
even  though  our  formulas  become  particularly  simple  when  the 
2nemn  are  small,  this  assumption  is  not  essential  to  our  procedure. 
We  also  define  6$^  as  the  difference  (in  cycles)  between  the  mean 
value  of  the  phase  of  ray  n and  the  mean  total-field  phase. 


6 = ♦(c  ) - ♦(c  ) , 

n n A A 


1 < n < 3 . 


(21b) 


Note  that  neither  the  A^  nor  the  $♦  are  required  to  be  small. 

n n 

The  last  assumption  needed  for  our  total-field  approximation 

is  that  AA/A  be  sufficiently  small  to  expand  the  total-field 

equations  in  this  quantity.  Again,  as  with  the  parameter  2ne.  , 

mn 

even  relatively  large  values  of  AA/A  near  0.5  give  useful 
approx imations. 
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Consider  first  the  double-path  cases.  From  Eqs.  (17) , (20)  , 
and  (21) , we  find 

(A  + AA)2  = A2  + A2  + 2A  A cos  [2n  (6*  - 6*  + e )]  , (22a) 

ro  n in  n m n mn 

where  m and  n are  the  indices  of  the  two  rays  involved.  We 

_ 2 

expand  Eq.  (22a) , neglecting  terms  involving  (AA/A)  and  those 
with  (2tremn)  to  powers  higher  than  the  first.  Using  Eqs.  (17) 
to  simplify  the  result,  we  obtain 

AA(t)/A  = 2nemn(t)  (An/A)  sin  (2x64^)  , (m  ? n)  . (22b) 


An  approximation  for  AL/(t)  then  follows  from  the  first  equality  in 
Eq.  (20c).  Rewriting  Eqs.  (17)  using  Eqs.  (20)  and  (21)  leads  to 
the  formula 

,n.A4  2iti  (A4>  +A$  )/2  i2ir  ( 64>  +e  /2) 

(A  + AA)e2lTlA$  = e m mn 


i2ir  (64>  -e  /2) 

+ Ane  n mn  ] . 


(23a) 


Expanding  Eq.  (23a)  and  simplifying  with  Eq.  (22b)  implies  the 
approximation 


A$(t)  = A*m(t)  - eran(t)  (An/A)cos  (2n6*J  , (m  ^ n)  . 


mn 


n 


n 


(23b) 


Thus,  amplitude  and  phase  changes  of  the  total  field  are  predicted 

A 

to  possess  the  same  time  behavior  as  A(t),  since  each  A4>  , and  hence 

m 

A 

each  £mn»  is  proportional  to  A(t). 

We  found  very  good  agreement  between  predictions  from  Eqs.  (22b) 
and  (23b)  and  numerically-computed  type  I double-path  results, 
described  in  Sec.  III.  The  (a,R)  = (11°,  7 km)  case,  with  m = 1 and 
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n = 2,  is  typical,  and  numerical  calculations  for  A(t)  and  4>(t) 

were  discussed  in  connection  with  the  solid  curves  on  Fig.  4. 

-4  -4 

The  mean  values  are  = 1.72  x 10  , A2  = 1.31  x 10  , 64>^  = 0*14 

_ _ 4 

cycles,  64>2  = 0.75  cycles,  and  A = 1.09  x 10  . As  an  example,  at 

A 

the  minimum  of  A (at  t = 9 h)  , the  values  of  A4>^  and  A4> 2 predicted 
from  Eqs.  (19a)  and  (19b)  are  -0.32  and  -0.25,  respectively,  and 
are  virtually  identical  to  the  numerically-calculated  values. 

With  these  numbers  in  Eq.  (23b),  we  obtain  A4>(9)  = -0.32  cycles, 
compared  to  the  numerically-computed  value  of  -0.31  cycles. 
Similarly,  using  Eqs.  (22b)  and  (20c),  we  find  that  AL(9)  = 3.8  dB, 
while  the  difference  between  the  numerically-computed  values  of 
20  log^QA(9)  and  20  log^g  £ is  almost  identical.  Theoretical 
predictions  for  phase  and  amplitude  are  shown  as  dashed  curves 
on  Fig.  4;  specifically,  the  amplitude  curve  is  20  log^A  + AL(t) 
using  Eqs.  (20c)  and  (22b),  while  the  phase  curve  is  A4>(t)  from 
Eq.  (23b)  plus  1.0  cycles,  to  correspond  to  the  way  the  solid 
curve  was  plotted.  Wherever  the  theoretical  curve  is  not  shown, 
it  is  indistinguishable  from  the  numerical  results.  As  can  be 
seen  from  Fig.  4,  the  maximum  difference  between  the  theoretical 

A 

and  numerical  values  occurs  at  the  maximum  of  A(i.e.t  = 3 h) , since 
the  numerical  results  are  slightly  non-symmetric  about  t = 6 h. 
Nonetheless,  the  relative  error  in  the  theoretical  predictions 
for  amplitude,  as  an  example,  is  less  than  12%  for  all  time,  and 
typically  the  error  is  less  than  this.  A similar  comparison  is 
shown  for  portions  of  the  (a,R)  = (12°,  8 km)  curves  in  Fig.  4. 

Even  closer  agreement  is  observed  here. 

For  type  I triple-path  cases,  formulas  similar  to  Eqs.  (22) 
and  (23)  may  be  developed.  The  analog  to  Eq.  (22a)  has  six  terms 
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on  the  right  side,  and  with  the  same  assumptions  on  (AA/A) 

2 

and  the  (2ire  ) , it  can  be  shown  that 
mn 

AA  = 2Tre23(t)  A3  sin  (2tt5#3)  + 2ne21(t)A1  sin  (2tt64>1)  . (24a) 

The  analog  to  Eq.  (23a)  is  expanded,  assuming  sufficiently  small 
- 2 2 2 

(AA/A)  and  ( 2tt ) (e  + e . ) /9,  to  obtain  after  simplification: 

mn  niK 

A4>  = A4>2  ( t)  - e21(t)  (A3/A)  cos  (2tt64>1) 

- e23(t)(A3/S)  cos  (2tt64>3)  . (24b) 

We  note  that  in  the  manner  they  are  written,  Eqs.  (24)  represent 
direct  extensions  of  Eqs.  (22b)  and  (23b) . 

The  approximations  Eqs.  (24)  produce  quite  accurate 
predictions  to  the  numerical  results  in  every  triple-path  case  we 
examined.  Typical  comparisons  are  shown  in  Fig.  6,  with  the 

same  conventions  as  Fig.  4.  In  particular,  for  (a,R)  = (8°,  8 km), 

-4  -4  _K 

the  mean  values  An  are  1.49  x 10  , 1.01  x io  , and  9.64  x 10 

for  n = 1,2,3;  A = 1.66  x io  and  64»n  are  0.11,  -0.02,  and  0.67 

cycles  for  n = 1,2,3.  The  theoretical  values  from  Eqs.  (19)  for 

A$n(3)  are  0.42,  0.30,  and  0.31  cycles  for  n = 1,2,3,  compared  to 

calculated  values  of  0.42,  0.30,  and  0.32  cycles,  respectively. 

The  predictions  from  Eqs.  (24)  and  (20c)  for  A4>(3)  and  20  log1QA+AL(3) 

are  0.38  cycles  and  -80.1  dB,  while  the  numerical  values  are  0.39 

cycles  and  -81.0  dB.  As  before,  agreement  at  other  times  is  even 

closer.  Also  shown  in  Fig.  6 are  comparisons  for  the  type  II  case 

of  (a , R)  = (10°,  10.5  km) . 

A particularly  striking  feature  of  these  approximations  is 
their  success  when  the  2iremn  are  not  particularly  small  with  respect 


to  one.  For  the  example  just  discussed,  the  maximum  values  of 
2ne^2  and  ^iTe^  are  0.76  and  0.67  radians.  In  iact,  the  derivation 
of  our  results  suggests  that  they  could  be  expected  to  produce  good 
approximations  so  long  as  the  2iremn  are  merely  less  than  one,  rather 
than  significantly  smaller  than  one,  with  the  important  provision 
that  AA/A  remain  small.  This  condition  is,  of  course,  violated  in 
type  III  multipath  behavior  and  in  type  II  cases  at  points  when  rays 
are  lost  or  gained.  We  also  remark  that  this  type  of  approximation 
to  the  multipath  field  takes  advantage  of  environmental  variations 
producing  small  relative  changes  in  the  phases  of  individual  ray 
arrivals,  not  on  the  specific  assumptions  of  the  environmental  or 
acoustic  model  in  this  paper.  Only  the  mean  phase  and  amplitude 
of  each  ray,  the  mean  total  field,  and  expressions  for  the 
individual  ray  phase  perturbations  are  required.  The  forms  of 
Eqs . (24)  suggest  possible  generalizations  to  problems  in  which 
the  total  field  consists  of  more  rays.  Such  approximations  may 
prove  useful  in  other  problems  where  conditions  for  their  validity 
are  satisfied. 

V.  SUMMARY 

In  this  paper  we  investigate  the  combined  effects  of  sound- 
speed  fluctuations  and  variable  bottom  depth  on  oceanic  sound 
transmission.  Off-shore  propagation  ranges  are  modeled  by 
assuming  a horizontal  ocean  surface,  a uniformly-sloping  bottom, 
a mean  bilinear  sound-speed  structure,  and  bottomed  source  and 
receiver  (above  and  below  the  SOFAR  axis,  respectively) . In  this 


I 
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paper  we  consider  propagation  in  the  down-slope  direction. 

Intervals  of  allowable  values  for  the  bottom-slope  angle  a and 
source-receiver  distance  R are  considered.  Sound-speed  fluctuations 
are  assumed  to  be  due  to  the  dominant  upper-ocean  temperature 
variations,  as  in  a previous  paper. ^ The  model  assumptions  guarantee 
that  rays  consist  of  planar  piecewise-circular  arcs.  For  our 
particular  assumptions,  it  is  found  that  for  each  pair  of  integers 
[K,L],  there  is  at  most  one  ray  having  K surface  and  L bottom 
reflections.  However,  the  set  of  SRBR  and  RBR  rays  that  does 
occur  changes  dramatically  depending  on  parameter  values.  An 
example  of  some  types  of  changes  is  given  by  Fig.  3,  which  shows 
regimes  of  existence  of  mean-state  rays  in  the  (a,R)  plane,  using 
fixed  values  of  other  parameters  and  a source-directionality 
condition.  Based  on  this  figure,  subsequent  investigations  are 
restricted  to  parameter  values  for  which  the  received  total 
field  has  at  most  three  contributing  rays.  It  is  emphasized  that 
other  model  assumptions  can  easily  enrich  the  possibilities  for 
occurrence  of  rays,  and  the  effects  of  factors  such  as  different 
mean-state  parameter  values,  different  source  and  receiver 
positions,  and  piecewise-linear  bottom  models  are  discussed 
briefly. 

With  expressions  for  travel  time,  spreading  loss,  and  boundary 
loss  and  phase  shifts  for  the  three  rays,  and  including  sound- 
speed  fluctuations,  the  total  field  is  computed  numerically  as 
a function  of  time.  All  single-path  cases  have  amplitude 
virtually  constant  in  time  and  show  a simple  proportionality 
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between  received  phase  and  sound-speed  fluctuations.  In  contrast, 
three  distinct  patterns  for  the  multipath  cases  are  categorized, 
each  of  which  may  have  significant  amplitude  variations.  We 
denote  by  type  I the  majority  of  time  evolutions  observed,  with 
roughly  sinusoidal  phase  and  amplitude  curves,  the  latter  showing 
peak-to-peak  variations  of  up  to  about  10  dB.  In  type  II 
behavior,  the  sound-speed  fluctuations  and  sloping  bottom  conspire 
to  produce  the  disappearance  and  later  reappearance  of  a ray  over 
a portion  of  the  twelve-hour  fluctuation  period,  leading  to  a 
jump  in  the  amplitude  and  phase  curves.  This  can  occur  for 
parameter  values  close  to  or  on  the  curves  in  Fig.  3.  A third 
type  of  amplitude  and  phase  characteristics  is  possible,  in 
which  the  time  behavior  differs  sharply  from  that  of  the 
fluctuations.  Type  III  behavior  arises  for  special  phase  and 
amplitude  relationships  in  the  mean  rays. 

A theoretical  explanation  is  developed  for  multipath  behavior 
of  type  I and  of  the  continuous  portions  of  type  II.  Time 
variations  of  individual  ray  amplitudes  can  be  shown  to  have 
negligible  effect  on  total-field  time  dependence.  Thus,  the 
first  step  is  to  relate  the  phase  variations  of  individual  rays 
to  the  environmental  fluctuations  and  the  variables  describing 
the  mean  ray  geometry.  Level  curves  of  maximum  phase 
perturbations  are  then  constructed  and  discussed.  The  second 
step  in  the  theory  is  to  relate  total-field  variations  to  the 
ray-phase  perturbations.  Formulas  for  the  deviations  to  the 
total-field  phase  and  amplitude  are  derived.  Then,  predictions 
from  the  formulas  are  compared  with  numerical  calculations.  As 
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illustrated  by  Figs.  4 and  6,  the  comparisons  show  quite  close 
agreement.  The  good  accuracy  of  the  theoretical  formulas  is 
discussed,  in  view  of  the  conditions  assumed  for  their 


derivation. 

In  future  work  the  authors  plan  to  extend  the  analysis  to 
include  other  combinations  of  environmental  effects  and 
acoustic  configurations,  including  depth-dependent  currents  and 
different  source-receiver  locations  and  ocean  boundary  models. 
For  any  such  problems,  characterizations  similar  to  the  two  in 
this  paper,  of  ray  occurrence  in  terms  of  parameters  and  of 
types  of  total-field  behavior,  would  be  needed.  Moreover,  it 
is  expected  that  the  theoretical  tool  developed  for 
approximating  total-field  variations  would  be  useful. 
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FIGURE  LEGENDS 


FIG.  1 


FIG.  2 


FIG.  3 


FIG.  4 


FIG.  5 


FIG.  6 


FIG.  7 


. Geometry  of  the  propagation  range  and  of  three  rays  of 
principal  interest. 

. Mean  bilinear  sound-speed  profile  and  pivoting  produced  by 
fluctuations. 

. Existence  regimes  of  various  rays  as  functions  of  range  R(km) 
and  bottom-slope  angle  a (degrees) . See  text  for  explanation 
of  numbered  curves. 

. Relative  total-field  amplitude  (dB)  and  phase  (cycles) 

versus  time  (h)  for  mean  double-path  cases:  (a,R)  = (11°,  7 km) 
and  (a,R)  = (12°,  8 km).  Solid  curves  computed  numerically, 
dashed  curves  from  theory  in  Sec.  IV. 

. Relative  total-field  amplitude  (dB)  and  phase  (cycles)  versus 
time  (h)  for  (a,R)  = (11°,  15  km). 

. Relative  total-field  amplitude  (dB)  and  phase  (cycles)  versus 
time  (h)  for  mean  triple-path  cases:  (a,R)  = (8°,  8 km)  and 
(a,R)  = (10°,  10.5  km).  Solid  curves  computed  numerically, 
dashed  curves  from  theory  in  Sec.  IV. 

. Level  curves  of  phase-perturbation  coefficient  of  sound-speed 

deviation  8 (cycles  m * sec)  in  (d-d  , R)  plane,  for  ray 
n r o 

(a)  n=l , (b)  n=2,  and  (c)  n=3.  See  text  for  explanation 
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period  in  the  upper  ocean.  An  examination  is  made  of  possible 
types  of  rays  for  down-slope  propagation  that  might  exist, 
depending  on  bottom-slope  angle  and  source-receiver  separation. 
The  total  acoustic  field  is  investigated  for  its  dependence  on 
these  parameters  and  time.  For  certa.infccnditions  when  up  to 
three  rays  comprise  the  mean  total  fieRl,  three  patterns  of  time 
evolution  are  described,  each  of  which  may  have  significant 
amplitude  variations.  Numerically-computed  examples  of  each 
type  are  presented.  The  linear  relationships  between  phase 
variations  of  individual  rays  and  the  sound-speed  fluctuations 
are  derived.  Then,  formulas  are  developed  to  explain  the  most 
frequent  behavior  of  the  relative  amplitude  and  phase  of  the 
multipath  total  field.  Predictions  from  the  formulas  show  very 
good  agreement  with  the  numerical  calculations. 


